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Abstract. Low-energy dipole excitations have been investigated theoretically in 
N=50 88 Sr and 90 Zr, several N=82 isotones and the Z = 50 Sn isotopes. For this 
purpose a method incorporating both HFB and multi-phonon QPM theory is applied. 
A concentration of one-phonon dipole strength located below the neutron emission 
threshold has been calculated in these nuclei. The analysis of the corresponding 
neutron and proton dipole transition densities allows to assign a genuine pattern to the 
low-energy excitations and making them distinct from the conventional GDR modes. 
Analyzing also the QRPA wave functions of the states we can identify these excitations 
as Pygmy Dipole Resonance (PDR) modes, recently studied also in Sn and N=82 
nuclei. The results for N = 50 are exploratory for an experimental project designed 
for the bremsstrahlung facility at the ELBE accelerator. 
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1. Introduction 

The progress in nuclear structure physics is closely connected with the extended 
experimental possibilities of the new facilities pp. In particular the experiments with rare 
isotope beams give us the opportunity to investigate nuclei far from stability. One of the 
most interesting results was the discovery of a new dipole mode at low-excitation energy 
in nuclei with high isospin asymmetry [HO [3]. Typically, one observes in nuclei with 
a neutron excess, iV > Z, a concentration of electric dipole strength of predominantly 
electric character at or close to the particle emission threshold. Since this bunching of 1~ 
states resembles spectral structures, otherwise known to indicate resonance phenomena, 
these states have been named Pygmy Dipole Resonance (PDR). However, only a tiny 
fraction, less than 1% of the total Thomas- Reiche-Kuhn energy weighted dipole sum 
rule strength is found in the PDR region. These studies have various astrophysical 
applications, e.g. the explosive nucleosynthesis and the neutron stars. 

Here, we present our investigations on the dipole excitations in many nuclei from 
N=50, 82 and Z=50 regions. For this purpose a method based on Hartree-Fock- 
Bogoljubov (HFB) description of the ground state is applied [6]. The excited states 
are calculated with the Quasiparticle-Phonon Model (QPM) jl]. 

2. The Model 

The model Hamiltonian [I]: 

H = H MF + Hf I+ Hf M + HZ (1) 

is built from the HFB term Hmf = H sp + H pair containing two parts: H sp describes 
the motion of protons and neutrons in a static, spherically-symmetric mean-field, taken 
as a Wood-Saxon (WS) potential. The parameters of the WS potential are derived 
from fully microscopic HFB calculations of the ground state [5l El [9], separately for 
every nucleus under consideration, which is different from the standard QPM scheme 
given in [3]; H pa i r accounts for the monopole pairing between isospin identical particles 
with coupling constants extracted from the data [TJ. The last three terms present the 
residual interaction H res = + H$ M + and refer to the multipole-multipole H V M 
and spin-multipole interactions H V g M of isoscalar and isovector type in the particle-hole 
and multipole pairing in the particle-particle channels. 

In the QPM the residual interaction is taken in a separable form 

^(r 1 ,r 2 ) = « A 3fi A (r 1 )» A (r 2 ) J (2) 

where 3? A (r) is a radial form factor , which is usually chosen as r A ; k x = (k X q, ft A i) 
are empirical isoscalar and isovector coupling constants, which are obtained by a fitting 
procedure [Hj. 

The nuclear excited states are constructed of Quasiparticle-Random-Phase- 
Approximation (QRPA) phonons, defined as a linear combination of two-quasiparticle 
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creation and annihilation operators as follows: 

Qt* = o E E (^^(jj'r) ~ tffA^VfT)) , (3) 
z ^ ii' 

where and are time-forward and time-backward^l operators, coupling proton 
and neutron two-quasiparticle creation or annihilation operators to a total angular 
momentum A with projection jj, by means of the Clebsch-Gordan coefficients C^, m , = 
{jmj'm'\Xfx). Correspondingly, 

At (jijoq) = \a + a + } = V C v ■ a + a + (4) 

^ raim2 

The QRPA phonon operators obey the equation of motion 

H,Q+\=E a Qi , (5) 

which solves the eigenvalue problem, giving the excitation energies E a and the 
wave functions of the excited states, defined by the time-forward and time-backward 
amplitudes ip^, and y^]')> respectively. 

The spacial structure of a nuclear excitation becomes accessible by analyzing the 
one-body transition densities 5p(f), which are the non-diagonal elements of the nuclear 
one-body density matrix. In the QRPA theory the one-phonon transition density is 
given by the coherent sum over two-quasiparticle transition densities entering in the 
wave function of a phonon by the relation: 

PA,(r) = E , (6) 

ji >h 

where the radial parts are given by the radial single particle wave functions and reduced 
matrix elements 

Pinq( r ) = R U r ) R nq(r)l(3i\^ X Y x \\j 2 ) , (7) 

A 



with A = V2A + 1. The BCS quasiparticle properties and QRPA state amplitudes are 
contained in 

<lf,, ''7.'/^ (V j2 +u j2 v h ) . (8) 
3. Application to PDR Excitations 



The calculated neutron and proton ground state densities are presented in FigJT] 
for N=50 and Z=50 nuclei. Of special importance for our investigations are the surface 
regions, where the formation of a skin takes place. For the N=50 nuclei, the neutron skin 
decreases from 88 Sr to 90 Zr, when the number of the protons increases. In the Z = 50 Sn 
isotopes we find, that for A>106 the neutron distributions begin to extend beyond the 
proton density and the effect continues to increase with the neutron excess, up to 132 Sn. 

X The time reversed operator is defined as = (— ) x ~ fJ "A\- fJ _ 
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Z= 50 




Radial distance [fm] Radial distance [fm] 



Figure 1. Ground state densities of N=50 and Z=50 isotopes used in the QPM 
calculations. 



Thus, these nuclei have a neutron skin. The situation reverses in 100_102 Sn, where a tiny 
proton skin appears at the nuclear surface. We find, that the properties of the ground 
states reflect directly off the low-energy dipole excitations. From QRPA calculations in 
N=50, N=82 [10] and 112_132 Sn [6j [9] nuclei a sequence of low-lying one-phonon dipole 
states of almost pure neutron structure, located below the particle threshold is obtained. 
The analysis of the dipole transition densities at E*<8 MeV in Z=50 (FigJUleft), N=82 
(Figj3]left) and at E*<9 MeV in N=50 (Figj2]left) reveal in-phase oscillation of protons 
and neutrons in the nuclear interior, while at the surface only neutrons contribute. 
These states we have identified with a neutron PDR. The states in the region E*= 
8-8.5 MeV in Z=50 and N=82 and E*= 9-9.5 MeV in N=50 nuclei carry a different 
signature, being compatible with the low-energy part of the GDR. At E*= 9-20 MeV a 
strong, isovector oscillation, corresponding to the excitation of the GDR is obtained. An 
interesting observation is the most exotic 100 Sn nucleus, where at E*=8.3 MeV a state 
with a proton structure is found. The analysis on dipole transition densities for different 
excitation energy regions in 100 Sn is presented in Figl5], illustrating the proton surface 
oscillations at E*<8.3 MeV. This mode could indicate a proton PDR. The dependence 
of the calculated total PDR strength on the mass number in N=50,82 and 10 °- 132 Sn is 
compared to the relative difference between the neutron and proton rms radii 

5r=yJ< >- v /<r2> (9) 

in the right hand side part of Fig.2, Fig.3 and Fig.4, respectively. In the case of 
N=50 and N=82 isotones we keep the neutron number fixed and change the proton 
number only. This affects the thickness of the neutron skin (see FigflJ left) as well 
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Figure 2. Dipole one-phonon transition densities in N=50 nuclei (left). The total 
PDR strength is compared to the nuclear skin thickness Sr, eqOH as a function of the 
mass number in N=50 nuclei (right). 



and respectively the total PDR strength (FigfS] right) and (Figj3] right) decreases with 
increasing proton number. The results obtained for 100 ~ 132 Sn nuclei, where the neutron 
number increases from N=50 to N=82 are in agreement with these obtained for N=50,82 
isotones considered above. Accordingly, the total PDR strength increases (FigJH right), 
when 5r increases and correspondingly the neutron or proton skin thicknesses increase 
(FiglH right). 

4. Conclusions 

In the isotones with N=50, N=82 and the Z=50 isotopes low-energy dipole states, 
identified with PDR were obtained. A close connection between the total PDR strengths 
and the neutron skin thickness defined by the relative difference of neutron and proton 
rms radii was found. These observations agree very well with our previous results for 
the Z = 50 isotopes and the iV = 82 isotones. In the most exotic nuclei 100 ~ 104 Sn 
lowest dipole states of almost pure proton structure are identified. They are related to 
oscillations of weakly bound protons, indicating a proton-driven PDR. The interesting 
point is, that these states are predicted to exist in heavy nuclei with N slightly larger 
or equal to Z. We suggest, that the effect is due to Coulomb repulsion, that pushes the 
weakly bound protons orbitals into the nuclear surface. The results for Sn isotopes and 
N=82 nuclei are in a good agreement with available data pUl [TT], [12] . 

At present, extended investigations on the fragmentation pattern of the low- 
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E*= 0-8.1 MeV E*= 8.1 -8.6 MeV E*> 8.1 MeV 




A 



Figure 3. Onc-phonon dipole transition densities in N=82 nuclei (left). The total 
PDR strength is compared to the nuclear skin thickness Sr, eqOH as a function of the 
mass number in N=82 nuclei (right). 




Figure 4. Dipole one-phonon transition densities in Z=50 nuclei (left). The total 
PDR strength is compared to the nuclear skin thickness 5r, eq02 as a function of the 
mass number in Z=50 nuclei (right). 
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Figure 5. QPM results for the one-phonon dipole transition densities in 100 Sn. 



energy dipole excitations are in progress. The QPM calculations will be performed in 
considerably larger configuration spaces and using microscopically derived interactions, 
thus enabling a detailed description of data on the dipole response of stable and exotic 
nuclei to be expected for the near future from ELBE and the experiments planned at 
GSI and for FAIR. 
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